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GBjECT 

It  Is  desired  to  study  the  orbital  stability  of  a 
space  vehicle  at  a  Trojsn  point  in  cis lunar  space*  i.e., 
where  the  earthy  moos  and  vehicle  form  an  equilateral 
triangle  in  the  orbital  plane  of  the  earth  and  220031.  If 
a  vehicle  originally  placed  at:  the  Trojan  point  is  indeed 
in  a  stable  orbit  relative  to  the  Trojan  point A  it  should 
then  be  determined  what  order  of  inaccuracies  of  injec- 
tion can  be  tolerated  to  still  establish  a  stable  orbit. 

This  assumes  that  there  is  an  operational  require- 
ment for  a  space  vehicle  of  this  sort  which  maintains  its 
relative  position  with  respect  to  the  earth  and  moon.  The 
requirement  might  be  for  a  space  navigation  or  cof&£unic&« 
tions  satellite,  ©r  it  might  be  for  a,  replenishment  station 
for  a  manned  deep  space  probe.  However,  the  specific 
advantages  of  a  Trojan  point  orbit  ware  not  investigated 
further. 


mmomosxa® 


One  solution  of  the  classical  problem  of  three  bodies 
Is  the  'equilateral  triangle  solution  discovered  by  Lagrange 
1©  1772 ,    re  the  three  bodies  are  each  equidistant  frega 
the  other  two  and*  in  the  absence  of  external  perturbations, 
maintain  this  relative  position.   Xn  the  restricted  three- 
body  problem  where  the  mass  of  the  third  body  is  small  enouv 
not  to  influence  the  motion  of  the  other  two  bodies,  the  two 
equilateral  triangle  points  in  the  plane  of  rotation  of  the 
two  saajor  bodies  are  called  Trojan  points.  The  Trojan 
asteroide,  the  first  of  which  was  discovered  in  1906*  have 
been  found  to  ramsin.  in  the  vicinity  of  the  equilateral 
triangle  point®  with  respect  to  Jupiter  m&  the  sun. 

The  term,  orbital  stability*  is  used  in  she  sense  that 

a  body   in     -jble  orbit  in  tl    rlnity  ©f  a  £roja&  point 
remains  indefirv    ■  in  this  vicinity,  i.e.,  there  is  an 
upper  bound  to  its  distance  fr*      -'.rojan  point  which 

. 

If  this  definition  is  accepted  lite         y,  we  sxust 
also  lisslt  tl  ^nitude  of  the  upi-  •tee  any 


clslunar  vehicle  w££fe  lass  thaa  esc&ps  velocity  will  remain 
In  clslunar  space,  barring  collisions  with  soother  body,  and 
will  have  an  upper  bound  to  Its  displacement  from  any  point 
in  the  earth-moon  coordinate  system,  h  cossvenient  restric- 
tion on  the  rigorous  definition  of  stability  la  not  immediately 
apparent.  However,  the  distinction  is  perhaps  overly  meticu~ 
lous.  We  are  concerned  with  a  vehicle  which  regains  in  the 
vicinity  of  a  Trojan  point  and,  due  to  initial  displacement 
from  the  Trojan  point,  initial  non~sero  relative  velocity 
with  respect  to  the  Trojan  point,  or  partubations  asternal 
to  the  system,  appears  to  oscillate  about  the  Trojan  point. 
If  it  continues  to  do  this,  it  is  stable;  if  it  diverges  from 
the  Trojan  point,  it  is  unstable. 

Historically,,  the  Trojan  points  have  been  considered  sta- 
ble for  "small  disturbances"  based  on  a  linearisation  of  the 
equations  of  motion  of  the  third  body  at  the  Trojan  point. 
However,  the  stability  of  a  body  in  the  actual  non*  linear  system 
has  not  been  conclusively  determined.  In  fact,  there  are 
arguments  indicating  that,  oven  granted  ezsact  knowledge  of  the 
physical  constants  &nd  unlimited  precision  and  speed  of  digital 
computation,  the  stability  can  still  not  be  established  with 
mathematical  rigor. 

The  existence  of  the  Trojam  asteroids  is  not  as  good  an 
argument  for  stability  as  it  seems  at  first  thought.  According 
to  the  solution  of  their  linearized  equations  of  motion*  the 
periods  of  their  two  modes  of  vibration  are  11.86  years  and 


467  years.  The  corresponding  pariods  in  eis lunar  space  are 
28.6  days  and  91.7  days.  Thus  the  54  years  available  for 
observation  of  the  Trojan  asteroids  corresponds,  in  one 
sense,  to  only  about  11  days  observation  of  a  cislu&ar 
Trojan  point  satellite. 

We  would  like  to  know  specifically  the  ranges  of  initial 
conditions  and  magnitudes  of  perturbations  for  which  a  Trojan 
point  vehicle  would  be  stable.  However,  the  study  described 
here  does  not  reveal  such  &  criterion.  Rather ,  it  seems  we 
oust  be  content  with  specifying  a  tiss  lisiit  during  which  the 
vehicle  will  not  exceed  a  specified  displacement  frost  the  Tro- 
jan point.  The  intent  of  this  investigation  was  to  provide 
answers  to  questions  of  the  sort,  "What  are  the  guidance 
accuracy  requirements  for  establishing  a  cislunar  Trojan  point 
space  station  to  remain  within  20,00J  miles  of  the  Trojaa  point 
for  at  least  10  years?"  Again,  Is  must  be  admitted  that  the 
answer  is  not  obtained  in  the  present  investigation.  However, 
the  directions  toward  which  such  answers  lie  may  be  indicated. 
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CHAPTER  2 


RESTRICTIONS  ON  THE  HtOBLBM 
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The  problem  Is  considered  to  be  that  of  the  restricted 
three-body  problem,  where  the  mess  of  the  third  body^  the 
vehicle,  is  negligible  compared  with  the  other  two  bodies. 
The  vehicle  will  then  have  np  effect  on  the  motion  of  the 
other  two  bodies .  This  is  obviously  valid  for  the  problem 
of  earth,  moon  end  a  vehicle  launched  from  the  earth* 

Another  restriction  is  the  elimination  of  all  pertur- 
bative  forces  external  to  the  earth-moon  system.  This  is 
equivalent  to  considering  the  barycenter  of  the  earth-moon 
system  as  an  inertia!  point,  and  neglecting  the  curvature 
and  velocity  variations  of  its  orbital  motion.  The  pertur- 
bative  effect  of  the  sun  is  perhaps  essential  in  the  solu- 
tion of  the  physical  problem.  However,  it  can  later  be 
added  to  the  differential  equations  as  a  forcing  function 
when  the  simpler  system  is  more  fully  understood.  Other 
perturbative  forces  such  as  that  due  to  Jupiter,  tidal 
forces,  earth  ob lateness  and  radiation  pressure  are  probably 
indeed  negligible. 

Also  the  orbit  of  the  moon  about  the  earth*  or,  mere 


( 


precisely*  about  the  baxycester  of  the  eareh-iaoom  system, 
is  considered  to  be  circular.  This  effect  can  also  be 
added  later  if  desired. 

Errors  due  to  these  sisaplifying  assumptions  3?111  be 
discussed  further  in  the  text.  Their  effect  is  also  con- 
sidered by  Bucheisi  (i). 
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CHAPTER  3 

THE  LXBRATXOif  fOIMB 


3.1  The  Geoaiatrv  of  a  Space  Vehicle  Attracted  by  Two  Bodies 

Zn  Fig.  3*1,  the  vehicle,,  m^,  is  attracted  by  the  two 
celestial  bodies,  su  and  a**  ^hieh  revolve  about  their  inter - 
tial  barycenter,  3*  In  elsiunar  space,  m^  and  m^  are  the 
earth  and  moon,  and,  for  convenience,  they  will  be  so  named. 
The  mass  of  the  vehicle  is  considered  small  enough  to  have 
negligible  effect  on  the  motion  of  the  earth  and  moon* 

To  simplify  the  equations,  all  distances  ere  normalised 
with  respect  to  the  earth-moon  distance.  Distances  from 
earth  and  moon  to  their  barycenter  are  then  e  and  \   respec* 

tively  'tikb&z® 

H        as* 


e  +  *  m   i 

The  gravitational  forces  on  the  vehicle  are  represented 
by  1%   and  1L  which  are  vector  specific  forces,  or  force  per 
unit  mass  of  the  vehicle* 

f,  -G»l  .eC^tej)/  »!  A  a  xO  ^i*^)     (3.1) 
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f    ^  „  G  frfr-k^)  /   «2   1  »  M^a^)   (3.2) 
r*        r2      *  2        r2 

Xt  is  convenient  to  choose  a  rotating  coordinate 
system  in  which  the  earth  and  moon  appear  stationary* 
The  equations  of  action  in  this  coordinate  system  must 
then  include  centripetal  and  Coriolis  accelerations 
resulting  from  the  coordinate  transformation.  Xt  is  also 
convenient  to  consider  the  centripetal  acceleration  of  the 
vehicle  due  to  coordinate  transformation  as  a  centrifugal 
specific  force ,  T ,  along  the  radius  vector  from  the  bary~ 
center • 

where  Q  is  the  angular  velocity  of  the  coordinate  system 
with  respect  to  inertial  space*  Then  from  Kepler's  third 
lav 

ft2«c  (m^) 


*inMvmmMmaaB*TWWM* 


(1)T  (3.3) 

f  «  G  (mfltBg)  r 

3 . 2  The  Existence  of  Librafcion  Points 

The  following  development,  due  to  Klemperer  and 
Benedikt,  (2)  shows  the  existence  of  five  points  in  the 
plane  of  the  moon's  orbit  at  which  a  vehicle  with  ssero 
relative  velocity  in  the  rotating  coordinate  system 
experiences  sere  resultant  force.  These  points  were 
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first  discovered  by  Lagrange  la  1772  as  particular 
solutions  of  the  three-body  problem.  Their  existence 
could  be  proved  with   less  effort  by  simple  substitu- 
tion of  tbelr  coordinates  into  the  equations  of  motion, 
but  the  following  development  also  show  that  these 
are  the  only  five  points  at  which  this  condition  occurs. 

Consider  the  vehicle  to  be  at  rest  in  the  rotating 
system*  With  sero  relative  velocity,  the  Coriolis 
acceleration  is  sero  and^  if  the  vehicle  is  in  equili- 
brium 

Consider  first  the  circuaf erential  equilibrium  of 
the  vehicle,  or  the  stmsaation  of  force  components  per- 
pendicular to  the  radius  vector  from  the  barycenter. 

-fx  sin  x  +  ^2  u%a¥'m  ° 

.     |  4   sin  x  *  — ™|-~*~  31^  Y*°  0     <3 A) 

-•4  sin  x  •  —9    &*&  & 
2 
One  locus  of  circumferential  equilibrium  is  the  straight 
line  through  m^  and  mg  where 

sin  x  •  sin  W    m  0 

Also.,  by  the  law  of  sines 
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£  i TTi  i«.i>Tii'-»i  ii  m i» mii  «-|  I  1 1  *  ■—  wiiw  iw'i 

rl  el 

Substituting  for  sin^and  sin  /^  is  Bq.  3.4 
*?  (€     sW )   ^         €   C?si»  <P  ) 

»m  iiiwijiMMMj/tiiaiii  T<wJfa»**ao      Suit  '»■  '■»■-  "Tidiw^iiiiwi  «»^»  Jp«u.>V:ifawi 

And  the?  other  locus  of  circumferential  equilibrium  is 
the  line 

Hesifc  consider  the  radial  equilibrium  of  the  vehicle* 
or  the  suasaation  of  force  ccssponests  parallel  to  the  radiue 
vector  frcss  the  harycenter. 

-f^  eos^  -  fg  coa^  +  f  »  g 

-  2£U^J — _    9om/T.  J£LJ£££L~  cos  <^+  rG  (s^tag)  «  o 

*1  r2 

cos  2f  -         e       cos  $  *$»  r  •  G 


1  2 


(3'5) 

Fln?3is!g  cos  X^n^  eosFoy  the  law  of  cosines 


e2  ©  t*2  *  r|  „  2*^  cos 
rl 


1     ""25     OT1 

^2    «.  r2.£,  y|  «.  2rr2  cos 


11 


2 
*2      2r~      ~2rr, 


coa^»r      +r2 


'2 

Substituting  for  cos  ^£  and  cosmic  Sq.  3.5 

Multiplying  by  -2^rJ  r|  to  clear  the  fractions 

?r2r2  +^rxr2  -/?e2  rj[  +  €r2r^  +  er^r|  -eA2**  -2r2r^r23    -  J 

(3.6) 
To  e-lindnate  r  frcss  the  equation,  r  is  found  in  tersse 

of  r-,  and  r9  by  the  law  of  cosines  . 


and  also 


r|    *  1  +  r2  -  2rx  cos    Jf  (3.8) 


Substituting  for1  coa^f  froxa  Eq.   3,7  into  Eq.  3.8 
r2  .  €2  +  r2  •  a  (X-fr^-r^)  •  e2-e  *  r^  (1-a)  +  eg2 

r2  »  «(l-a)  *  r|  (l-«)  +'  *r|  ■  «#  +  *r|  +  er2 

2 

Substituting  for  r     in  Eq.   3.6  and  combining  terms  the  lecua 

for  radial  equilibrium  is  obtained 

2       5  3^2         35       q  3  3        351+6       32       1+623 
-T-  *1*2  +  7F    rlr2  •  2rlr2  "  *X~r2~  X"  *lr2  -  T  rlr2  + 

3         3 

+r^  +  r2  «  0 

The  equilibrium  loci  are  shown  in  Fig.  3.2.  The  five 
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intersections  of  the  t&o  sets  of  curves  are  the  only  points 
at  which  a  body  would  remain  at  equilibrium.  These  five 
points  are  statically  stable  in  the  rotating  coordinate 
system,  but  the  development  yields  no  dynamic  information. 
To  specify  the  behavior  of  a  vehicle  in  the  vicinity  of  a 
Vibration  point  with  deviations  in  position  and  non-eero 
relative  velocity,  the  equations  of  motion  of  the  vehicle 
must  be  investigated* 
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Fig,  3,2.  Location  of  the  Libration  Points 
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CHAPTER  4 


lint  M»  'i  I  ■■■! 


the  tomtichs  of  msim 


4.1  Foxaaiil^fcioii  o£  the  Equations  of  Motion 

In  an  inertia!  coordinate  systsss.,  the  ©otlen  of  the 
vehicle  can  he  described  by  equating  its  acceleration  to 
its  gravitational  specific  force.  Since  it  is  desired  to 
describe  the  csotion  in  the  coordinate  set  vanish  rotates  with 
the  mxm,   the  additional  apparent  accelerations  which  appear 
in  the  coordinate  transformation  must  be  included.  The  total 
acceleration  consists  of  the  sum  of  the  following  terse®  s 

1*  The  acceleration  of  the  coordinate  system  in 
inertial  space.  IMs  is  considered  ser©  in  the  restricted 
three-body  problem.  In  the  physical  problem  in  cislunar  space, 

it  corresponds  to  the  acceleration  of  the  earth~moen  system  in 
its  orbit  around  the  sun,  including  the  perturbation  effects 
of  other  planets.  Its  only  significant  component  is  the  cen- 
tripetal acceleration  of  the  earth  which  seems  to  be  the  ssost 
appreciable  effect  neglected  in  confining  the  investigation  tc 
the  restricted  three-body  problem*  If  later  it  appears  desirable 
to  consider  this  term*  it  can  be  included  as  a  periodic  forcing 
function  in  the  equations  of  motion,  representing  the  differ* 
between  the  specific  gravitational  force  of  the  sun  on  the  - 

XL 


vehicle  and  that  force  on  the  barycenter  of  the  earth-moon 
system, 

2.  The  tangential  secelera&ic&«  This  is  also 
considered  to  be  ssro  since  the  moon  is  assumed  to  be  in  a 
uniform  circular  orbit  about  the  earth*  If  the  effect  of 
the  eccentricity  of  the  moon's  orbit  were  to  be  considered . 
this  term  would  have  to  be  included, 

3.  The  centripetal  acceleration*  This  is  treated 
as  a  centrifugal  specific  force  and  grouped  with  the  gravita- 
tional forces. 

4.  The  relative*  acceleration.  Xn  the  {%,  y)  rota- 
ting coordinate  system.,  its  components  are  k  and  y"  which  ars 
included  in  the  equations  of  motion* 

5.  The  Coriolis  acceleration.  Its  magnitude  is 
&    *$»  y  ,  and  its  direction  Is  defined  by  the  vector 

cross  product  of  the  system  angular  velocity  and  the  relative 
velocity  of  the  vehicle*  Its  x  and  y  coacponents  are  -2£ty  and 
+2S&  respectively. 
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Fig.  4ol.  Representation  of  the  Rotating  Coordinate  Set 


As  shma  in  Fig.  4.1,  the  origin  of  the  coordinate 
system  is  chosen  at  en  equilateral  triangle  point,  or 
Trojan  point,  since  t?e  are  pri&sarily  concerned  with  the 
dynamics  of  a  vehicle  in  the  vicinity  of  this  point. 

The  at  and  y  ccsaponenfcs  of  the  vector  equation  ©£ 
motion  are  then 

iy   <£y   cy 

A 

Substituting  frcns  equations  3-1,  3.2,  and  3.3  and  resolving 

x  -  £}£  .^(ai*^)  cos  y  +  €G  (y^ } cos  3  *  Gdaj+n^)  r  cos  § 

rl  2 

y  +  20&  -  ^G^l^>  sin  y +  €f<8alHa2>  sin  ^  +  6(Bu4flu)  r  sin  * 

1 
To  further  simplify  the  equations,  time  will  be  normalised 
with  respect  to  the  period  of  rotation  of  the  earth-moon  system 
so  that 

t2  „  OCy^)  t*2  w2iere  fc.  l3  now  real  time 
IF" 


2  2 

dt,z     A  *  at'* 
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with  the  dot  denoting  differentiation  with  respect  t$ 
dimensionless  time* 

In  dlmsnsionless  time,  the  equations  of  motion 
become 

*  *  /?  <•  si  dh 

x  -  2  y    »  -  -£±-  cosy+— -j  GO&F+  r  cos   /   (4#1) 
-»2  r 

r  2 

y    4-  2  x  »  -  — j-    sin  X  -  —2^  sin^-fr  r  sin  y 
rl  r2  <**2) 

Ref<?rring  to  Fig.  4.1,  the  quantities  rp  r,„   end  the 
trigonometric  functions  of y,  P  9  and  r  may  be  expressed 
in  terms  of  x  and  y  as  follows: 

r?  »  (y  +  ¥1)    +  (s  +  1)  •  y2  4-Sy  ♦<2  +  x+l 
1  .2.-2  .... 

r2«(y+    ~)    +  <--«)     •  y2  +^3y  ♦  x2  -  x  +  1 


r  cos    f«  x  >  /}  -    —  ■»  X  4*    — -c 

2  -     2 

r  sin  f  «    y+    - 

2 

x  +  %  or  'j  "% 

cos  Y  »    ~jt—    •  <*"  +  %)   (y    *  *3y  4-  vT  +  x  +  1) 

cos  /?  «.    JLZJS  -  -  (x  -  %)   (y2  4-  ^:ly  +  x2-  x  +  1) 
r2 

y+^l  fik       2     f         2  -% 
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=~  «=  (y  +  t,-)  (y  +  3y  ■ 

2 


in  ^  - 


Substituting  in  equations  4.1  and  4.2  the  equations 
of  motion  become 


x-2£  «  sefZg-e- 


(»*-?i)(y2+  3y+3i2+»H)/^2„e(2c-y  (y2-fjjrfx2-x+l)  '2 


(4.3; 

y"+  2x  «  y+  -  -  (yf  ~)(y2+  .jy4-x2+::+l) /^- :(y+  ~)  (y2+3y+x2-js*l)  /2 
2       2  2 

(*.4) 

4.2  Integration  of  the  Equations  of  notion 

The  only  l:novm  integral  of  those  equations  was  discovered 
by  Jacobi.  It  has  since  been  proved  that,  in  this  coordinate 
system,  no  other  algebraic  integrals  exist. 

Let  u  =  -cx+^(y2+jy4.::2;:+l)+^{y2+3y^::2+x->l) 


•f^y^^y+^-xfl) " 


(4.5) 

Then  —  «=  -efeeMs-zK^)  (y2+3yf-x2*x+l)  /2 

^2:  -3 

"€ (x~ij)  <y2+»6yfK2-:ftH)   2 


(4.6) 

3is  «  y  +\3  ..  7|  (yf  .2)  (y24Ty+x2+x*l)  /2 
ay      2       2 

-G(yf 

2 


'A.   Of        9     ~3/0 

„e(yf  ^)(y*+  3y+x  -:frH)  /z 
2 

(4.7) 

and  the  equations  of  motion  are 

S  -2y«  ^  (4.3) 

dx 

19 


y+2*o2s  {4.9} 

Multiplying  equations  4*3  and  4,9  hy  2x  and  2y  ^respectively 

2jrx-2yxa2a^  (4.10) 

dx 

2"yy  *  2yi~  2*  ^  (4.11) 

dy 

Adding  equations  4*10  m&  4.11 

2xs  *  2  y  y  *  3x  ^  +  2y  — 

&s  -    dy 

2x&+2}&«»2&    fe  *  2  .§£   ^ 
dt  ,dt  dx    dt  -       dy    dt 


c  •  -2u  +  i2+y2 

The  Jacob!  integral  is 

C  -  2€x  -  (y^^y^^lJ-a^Cy^^ss2^-!) 

-2«(y^/3j<«2-«*'lJ*  *  **  +  J*  •  constant 

(4.12) 
4*3  Vectorial  j£epre@antati@n  of  the  Equations  of  Ifotioa 

It  should  be  noted  bare  that  the  subsequent 
equations  apply  to  only  one  of  the  two  Trojan  points*  For 
the  other  Trojan  point*  a  syaaaetry  about  the  x-axis  applies, 
but  the  equations  and  directions  may  not  be  applied  literally 
without  consideration  of  this  syasaetry 

If  the  motion  of  the  vehicle  is  represented  in  the 
cocsplex  plane  where 
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s  «  s  4-  iy 


equations  4,3  and  4*4  can  be  written  as 


35  «  2  y  w  £. 


K 


(4*13) 


«i 


iy  4-  21s  •  if 
Adding  equations  4*13  and  4,14 


<4a4) 


s  *  Ztz  ®  f„  4-  if. 


(4*15) 


At  first  thought,  it  would  seem  that*  sor  a 
Trojan  point  orbit,  the  resultant  of  £x  and  fy  should  be 
directed  toward  the  origin,  Zhis  is  not  the  esse;  is  will 
be  seen  later  that,  in  the  vicinity  of  the  Trojan  point,  i% 
and  f  to/©  the  seese  sign  as  x  and  y  respectively.  Rather, 
it  is  the  Coriolie  acceleration  which  provides  an  average 
relative  acceleration,  z]     toward  the  origin  if  the  orbit  is 
to  be  a&able.  The  vectorial  .representation  of  Eq*  4* IS  is 
shown  in  Fig.  4,2*  This  indicates  that  the  orbit  must  have 
a  gross  motion  clockwise  about  the  origin  for  stability* 


Vtw«-»»TCT»l  MHWIIItUOMWaC**"    ""  "  "■■'"  "  '"■'  ■  ■  3f3*       ^ 


.  4.2.  ¥e 


Representation  of  the  Equations  of  Motion. 
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4 « 4  Linear  isat ion  of  the  Equations  of  Motion 

As  a  first  order  approximation  of  the  behavior  of  a 
vehicle  in  the  vicinity  of  a  libratlon  point ,  the  equations 
of  motion  are  linearised  and  the  solution  of  the  linearized 
equations  is  examined.  The  procedure  used  by  Moultcn  (3) 
is  essentially  to  linearize  the  right  hand  sides  of  equations 
4.1  and  4.2  by  Taylor  series  and  solve  the  resulting  linearised 
equations.  Whittaker  (4)  uses  a  different  approach,  working 
with  the  equations  of  motion  in  the  form  of  four  first  order 
differential  equations,  with  the  sacse  results. 

The  first  approach  will  be  extended  here  to  obtain  a 
more  complete  linear  model  of  the  system,  including  an  integral 
of  the  equations  equivalent  to  the  Jacob  tan*  This  will  later 
be  used  for  a  more  complete  comparison  of  the  actual  and 
linearized  systems. 

Linearizing  the  differential  equations  is  equivalent  to 
retaining  second  order  terms  in  the  Taylor  series  repre- 
sentation of  the  Jacobian  function*  Eq.  4-5,  Letting  primed 
symbols  represent  the  quantities  of  the  linearized  system 
Eq.  4.5  can  be  written  as 

u'  (x,y)  2  u  (0,0)  +3S  (0,0)xe.|g  (0,0)  y+%  £-8  <0,0)«*  + 

(4.16) 

-3?a  <0,0)Ky«*  &4   (0,0)y2 
dxdy         $y* 

Substituting  x  »  y  »  0  in  equations  4.5,  4*6  and  4.7 
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U  (0,0)  ■  \  -&3+   €  *  — 

2 

— §~  (0,0)89  -€  ♦  %  -  %^4»  %  £  =  0 

-£H  (0,0)  «  1  -^  -  €  -  •  0 
dy        2  n  1  2 

Different lacing  equations  4.6  and  4*7 

iLH-*»  l+3^(x4>5)2(y2+Sy4«2+2i?fl)   ^2  v7(y2*€y«-x2+x+l)     '* 
dx2 

+3s(x-%)2(y24-^3y^x2-^l)'"54  -€   (y^^yfrx2-^*!)    '2 

(4.17) 

-^«  l*3^(y*~)    <y2+^yl«2+^l)>?/(yVv53^24flBfl)"^ 
dy  2  , 

+3«(yfl5)2(y2+V53rte2  -»rt)T2€  {^^^.M)"3/i 

(4. 18) 

i-S  »  3d{xt-k)    (yff)  (y2*{iy*^2^lf /^  *  3G(x-%)(y^|)(y2^5«2-»}-ir2 
Sxdy  "  ^ 

(4.19) 
Substituting  x*  y  ■  0  in  equations  4«X7,  4.18  and  4.19 
2 

—. *  (<>*o)  «i4-~a-^~c-€B5~ 

fcx^  4       4      4 

2 

£-H  (o,o)  •  1  +  ^-^  +  2  €.€«  2 

dy^  4      4      4 

dxdy    '       4       4  4 

Substituting  these  quantities  in  equation  4. 16, the  Jacohian 

function  for  the  linearised  model  is 


23 


2    8         4  8 

Differentiating,  the  linearised  equations  of  tsotioa  see 
obtained 

i"  «  2  y  •   «&'  «  I  «  *  3*2  (i  -  2«)  y 
te    4   -   4 

(4.21) 

y*  +  2  i  •  ^'  -  -21^  (1-2S)  x  +  £  y 
&y     4         .4 

(4.22) 

end  the  integral  of  the  ©quatioias,  ©quiv&laofc  £0  the  J&cobtaa* 
18 

C»  »  x  2  +  y  2  -  3  «  ^  x2  -  -2^1  (X~2«)  xy  -  ^  y2 

4       2  & 


<4c23) 
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V 


CHAPTER  3 

•wcauBMuncaw 

ire  sazmtm  of  xbs  ubbarzzep  gocagioRS  or  motion 

5.1  Tha  Stability  Crlteglcro 

Xa  differential  operator  form,  equaticas  4*21 
and  4*22  can  be  wslztm  as 


4    *   k-     4      ' 

The  ehar&e&eri&tlc  equatiea  o£  this  ©ystesa  id 

<P2  -  ~)  <P2  -  -)  *  |2p  *  A-^  <1«W  2p  -  &■£  <X»2«)J»  0 

4       4    *—     4     — *       4 

p4  «§*  p2  +    ^2  €  (1«6)  «  0 

4 
the  roots  o£  the  eharacteri@sle  eqtaatioa  are 


If  tha  second  teres  is  seal*  the  characteristic  ec.uatioa 
yields  two  pairs  of  pure  imaginary  soots,  aad  the  solution 

coaeia&s  o£  pure  ate^soido*  Itos^er.  1£  She  eecosd  teas  is 
imaginary,  the  roof ©  have  poalt&re*  seal  parts,  a^d  the  sns- 
piitude  of  the  solution  increases  ©aponeafclally  with  titeae* 
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The  stability  cxttmtlmx  is  then 


OT 
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**»**>  •    %  -ii  -    ^  -  0.03852 


(I)        m    L-S-.)    m  0.04006  « 

max 


1. 

24.96 


'1 


In  the  Q&£th"W0Q&  system*  ™  is  approximately    ^srv*  ssut 


the  stability  e&i&erios  is  met. 

S*2  The  gtraiffat-Liae  Libration  Folate 

SfoultOQ  wg(sa  a  similar  de^alopgasat  to  show  that  the 
straight- Usae  li&ratioa  points  are  unstable.  Aaother  presen- 
tation of  this  approach  is  mede  by  Bueheim  (5). 

If  the  equations  of  motion  are  linearised  about  &  straight* 
line  libration  point*  £hs  charaeteristia  eisuatioa  of  the  result- 
lag  system  yields  tsro  imaginary  root©  ©ad  two  real  roog©5 
one  of  <a&icb  Is  always  positive*  A  part  of  the  solution 
tmst  thea  increase  exponentially  with  time*  The  system 
is  thus  eosasidoffed  unstable!  a  small  disturbance  ia  displacement 
from  the  libr&tion  poijit  will  increase  with  time* 
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5>3    The  Ygpjaa  Point:  Solufcloa 

»»i»iw>— laawniin  iwOHIwiimn ■»■»«»— — m—  k.iumi— »w>* 

nu 
If  fcfce  as$0s  &M&io    «~-    ie  stroll  eaoisgfa*  nit  if  ia  die 

■■'! 
ear£h««soeffi  ©y^fceaa,  fc&a  linearised  equeeioa®  of  moctoa  yield 

a  soiufcioa  of  ste  fosm 

X  •  A^  aim  p^£  «§*  B^  Goe  p^fc  *  G^  sin  pgfc  +  D^  cos  p£fc 

(5.1) 

y  *  Ag  sis  p%&  ♦  B2  000  p^fc  *  c£  s^°  P2fc  *  D2  cos  ?2fc 

(5.2) 
vbere 


aod  fouse  of  the  ©£gfe£  ©osstaofes  of  £&eegff«&io&  are  &ndepeade»&< 
X8w  relAfciosiflfeipa  seoog  Che  cowztmafc*  essa  be  defcemfcaad  by  aub« 
HfcifcutiUag  She  solu&ion  bask  issto  fcisa  llaeesr  differencial  equ&~ 

tioas*    Solving  for  efccs  eongtgf&s  of  Eq«  5*2  fra  fcems  of  She 
coaatante  of  Eq«  5*1 


■a 


11  Tin  ■  1 1  iniii—iTiii  «*"■  V«/««5J 
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%  "  ^  f$   (1-26)   B 


i 


pf+1  <5-4> 


Jf3  (l-2€)  cx  -  2?^ 
?2  +S 


Cl  '  ~  sft  (1.2t)  0X  (5.6) 


9-9  as 

P2+5 

The  constants  A,,  Bp  Cj,  Dj  are  than  determined  by 

the  initial  conditions  of  displacement  and  velocity,  x  (0), 
„  *  C©)  and  £  <0) . 

■  ■  ■■■'■        -      ";  ■  '  •  ■  '  ■■'■;-'-:  Q4';K J-- .'!■  ■  -■'■   &$■■*$  ■   '  '':c-  ":'\, 

to  better  understand  the  nature  of  the  solution  that 
lias  been  obtained  from  the  linearised  equations ,  it  would 
be  desirable  t©  deteraine  how  the  solution  is  divided  bei 
tbe  two  modes  of  vibration  as  a  xuactioa ' o£  initial  conditions. 
The  constant®  of  equations  5.1  end  5.2  can  be  determined"  as 
a  function  of  initial  conditions ,  in  terms  of  the  natural 
frequencies  p.,  and  p*«  This  is  doate  1b  Appendix  A  using 
numerical  const     £®r  the  earth-saoon  system,  but  the 
physical  process  b;  axclted  to  differ'; 

proportion©  o£  the  two-  asodee  of  vibration  is  not  readily 
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To  obtain  further  insight  along  this  line,  we  determine 

what  specification  of  Initial  conditions  will  excite  only 
one  G&ode  of  vibration,  say  the  p*  mode*     For  simplicity  we 
arbitrarily  specify  the  Initial  displacement  as 

x  (0)  «  0  (5.7) 

y  (0)  m  o  (5.8) 

We  are,  of  course,  concerned  with  much  smaller  displacements.. 
but  the  results  can  be  scaled  do^n.  If  we  specify  that  the 
initial  velocity  is  chosen  to  eliminate  the  p^  zaode,  then 
equations  5.1  and  5.2  become 


x  »  A*  sin  p^t  +  B,  cos 
y  •  &2   8*n  Pxfc  *  B2  cos  pl£ 
Substituting  initial  conditions 

X  (0)  to  R^  ra  0 

y  (0)   »  B2  •  1 
x  (0)  •  pjAj 
y  <0)  «  Pl^2 
Also  from  £q.  5.3 

Ag  »."  Jff  (l-2«)  Ax 
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\*)  ■ 


A+i 


Sine®  3*  *?  1       than  A,  »     1_  •  2  +  9,  ) 

2p       l      4 

and  Aj  •    -  -L      -2/3  0-2«* 
Thus  the  initial  velocities  ssust  be  specified  as 


*  (0)  «  h  <P?  -2-    -  )  (5.9) 

1        6 


3   — 

J73  (1-2*)  (5.L>) 


And  tho  subsequent  isoticm  of  the  particle  is  specified  as 
z  »    pl  +  4         ©in  ptt  (5.11) 

3 


y  e» 

*.,  I?,  CkZSl     ain  Plt     4-  coa  p,t       (5.12) 
Jpj 


y  CT  ?1C      *    CO*    P^ 


Solving  for  sin  pjt  and  cos  p^fc  in  terms  of  %  and  y,  we 
can  fehsn  eliminate  time  by   the  equation 

[stn2  P^t  I(x,y)  *  cos2  pjt  (x,y)  ■  .'. 

After  siBEplification,  the  equation  of  the  "orb:f.t"  about 
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the  Trojan  point  is  found  to  be 


fi+sSjkiEl      xy  A"   (l-2€)     »4pf]    2 

1      4  ll6<P?+*)       J 

4 

(5.13) 


This  is  the  equation  of  an  ellipse  with  its  center  at  the 
origin.  It  shoulc  be  noted  that  this  is  not  a  Kepler ian 
elliptical  orbit.  Angular  momentum  is  not  conserved  in 
this  coordinate  tyatem,  and  there  is  no  central  force 
field  analogous  f.o  that  o£  a  Kepler  ian  orbit  in  inert  ial 
space. 

Although  th;<  motion  of  the  vehicle  about  the  Trojan 
point  described  La  the  relative  coordinate  system  may  not 
be  correctly  terted  an  orbit,  the  term  will  be  retained 
for  convenience. 

If  we  let  v  be  the  magnitude  of  the  velocity  of  the 
vehicle  in  the  relative  coordinate  system,  i.e. 


v2  »  &2  -5-  f2 


Then  this  velocity  can  be  found  at  any  point  on  the  orbit 
from  Eq.  4.19 

v2  «*  0?  +  3  *  ~  x2  +  -2  f5  (l-2e)  xy  +  2  y2 
4      2  4 
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C  is  a  constant  which  can  be  calculated  from  initial  con- 
ditions 

C*  -  v2  (0)  -  3  -  ~ 

4 

Fro©  equations  5.9  and  SdO 

v2  (0)  en  %  (p|  +  |)     *  27  <X»2<02 

- 

*nr*  at  any  point  on  the  orbit 


v2  -  k   (pf  +  9)2  *  g  (i-2«)  -J  *  |  k2  +  |fj  (Me)  ^ 

♦-y2    (5.14) 


Considering  an  arbitrary  point  on  the  orbit  to  be  the 
initial  point,  the  initial  velocity  necessary  to  excite  only 
the  p«  mode  is  determined.  Its  magnitude  is  obtained  from 
Eq.  5.14,  and  its  direction  is  specified  by  the  tangent  to 
the  orbit  computed  from  Eq.  5.13. 

If  we  now  divide  all  displacement  and  velocity  varia- 
bles in  equations  5.7  through  5.14  hy  an  arbitrary  scaling 

factor,  \>ie  obtain  a  new  orbit  with  the  same  center ,  shape 
and  orientation,  which  is  also  a  solution  of  the  linearized 
equations  with  initial  conditions  specified  to  eliminate 
the  p2  mode.  For  a  continuous  variation  of  scaling  factor, 


32 


an  Infinite  f 

gaily 

of 

ellipses 

25(G)  »  )  0 

f(Q)  •      k 

k(Q)  «*      - 
2 

<pf  * 

9\ 
4 

£(0)   m  2k\$ 
8 

(1~2«) 

2Pl 

♦  2) 

-   4 

sin 

plfc 

covering  the  ss-y  pl&ng* 

(5,15) 

(5,16) 

(5.17) 

(s.ia) 


y  »  Sk  fe  (1-2C)  sin  p,  *  k  cos  p,t 


^.JjfLfldBL^    *y* 
2<pf  4.  J 


,2 


2 


2       9n 

7(1~2«)  *  4p£ 

6(pf*    S) 

A  „     4 


2 


(5. 19  J 

(5.20) 

(5.21) 


ciPi  + 1)  +  aa  (M€)*.  S4V  + 


V      «* 


64 


4        4 
£  f3(l-2«)xyt  -  72 


\a)  f.it&J 


divan  say  point  in  the  x,  y  plane  as  the  initial  point* 
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we  can  calculate  k  hy   substitution  in  Eq.  5.21.  We  can 
then  determine  the  required  initial  velocity  vector  to 
excite  the  p^  mode  only  from  equations  5.17  and  5.18  and 
the  resulting  orbit  is  specified  by  equations  5.19  through 
5.22. 

If  we  wi@h  to  eliminate  the  p-j  mode  instead  of  the 
?2  mode,  the  procedure  is  the  sase.,  and  the  resulting 
equations  are  identical  with  equations  5.7  through  5.23 
with  the  substitution  of  p«  for  p^. 

5*5  Determination  o£_tha _  Orbit  from  Specified  Initial 

Conditions 


We  now  have  the  mathematical  equipment  for  answering 
the  question,  "Oiven  an  arbitrary  initial  position,  what 
vector  initial  velocity  is  required  to  excite  only  mode 

pl  or  p2?*  We  ^°  eot  ^ave  GS  conveB^ent  a  physical  repre- 
sentation for  answering  the  question  "Given  an  arbitrary 
Initial  position  and  an  arbitrary  initial  velocity,  what 
does  the  orbit  look  like?"  The  temptation  is  to  substi-* 
tute  numbers  into  the  equations  of  Appendix  A.  However, 
we  can  do  it  semi-graphically  hy  a  two  stage  trial  and  error 
method,  and  perhaps  gain  a  better  physical  insight  of  the 
problem. 

Leg  ilis  equations  of  the  orbit  be  written  as 


34 


x  -  *l  +  x2  (5*23) 

*  "  yl  +  y2  (5.24) 

where  x^  -  A^  sin  pjt  +  B-^  cos  p,t 
vl  •  ^2  sin  pjt  +  B2  cos  p2t 
x2  •  C1  sin  p2t  +  Dx  cos  p2t 
y2  »  C2  sin  p2t  ♦  D2  cos  P2t 

To  determine  the  orbit  resulting  from  specified  initial 
conditions,  pick  a  point  (jt^O),  yi(°)J  as  a  trial*  This 
determines  an  ellipse  of  the  family  of  Eq.  5.21  and  specifies 
the  required  x^(0)  and  $^(0).  It  also  specifies  the  point 
I  *2^'  y2(0)J  by  equations  5.23  and  5.24.  In  turn,  this 
point  also  specifies  an  ellipse  of  the  family  of  Eq.  5.21 
with  p2  substituted  for  pj,,  and  gives  the  required  x2(c)  and 
^2(0).  If  the  trial  is  correct,  the  vector  addition  of  v«(Q) 
and  v2(o)  gives  the  specified  initial  condition  v(0),  and 
the  orbit  can  be  constructed  by  the  addition  of  a  vector 
generating  the  p«  ellipse  at  a  relative  angular  velocity  of 

El 

Pi  ' 

This  method  is  illustrated  in  Fig.  5-1  which  is  not 
drawn  to  scale.  It  is  included  only  for  possible  clarifica- 
tion of  the  above  explanation. 
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5.0  Desirable  Ranges  of  Initial  Conditions 

From  t^e  preceding  analysio,  it  appears  that,  given  an 
unavoidable  displacement  error,  we  will  obtain  the  "smoothest" 
orbit  if  the  initial  velocity  has  a  magnitude  in  a  range  such 
that  the  Jacobia»\ constant  of  the  vehicle  is  approximately  the 
same'  as  that  of  a  body  at  rest  at  the  Trojan  point   It  can 
be  seen  by  substitution  in  Eq.  4 .19  that   a  body  vibrating  in 
the  p,  mode  only  has  a  relative  energy  higher  than  that  of  a 
body  at  rest  at  the -Trojan  point,  while  a  body  vibrating  in 
the  po  mode  [only   has  a  relative  energy  lower  than  that  of  a 
body  at  rest  at  the  Trojan  point'. 

More  important,  we  wish  the  initial  velocity  to  be  in  a 
direction  compatible  with  a  clockwise  orbit  about  the  Trojan 
point  (using  the  notation  of  Fig.  4.1), 

The  asethod  of  the  preceding  section  shows  that,  for  ini- 
tial velocities  in  the  "correct"  direction,  we  tend  to  get 
contributions  from  the  two  modes  which  are  additive,  and  the 
outer  bound  of  the  resulting  orbital  motion  is  not  much  different 
from  the  initial  displacement.  On  the  other  hand,  if  the  ini- 
tial velocity  is  in  the  opposite  direction,  we  tend  to  get 
large  contributions  of  the  two  modes  which  subtract  at  the  ini- 
tial point,  but  &dd   at  the  outer  bound  of  the  orbit. 

This  verifies  the  impression  we  ge£  from  section  4.3. 
Since  the  gravitational  and  centrifugal  force  field  is,  in  general, 
directed  away  from  the  Trojan  point,  the  Cor io lis  acceleration 
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taust  provide  the  net  acceleration  toward  the  Trojan  point 
for  a  stable  orbit.  !7or  this  to  occur,  the  gross  motion  of 
the  vehicle  mast  be  clockwise  as  can  be  seen  in  Fig.  4-2. 

Fig.  5-2,  whicl  is  also  not  drawn  to  scale,  is  a  sketch 
of  two  attempted  in;,  act  ions  of  a  vehicle  into  a  Trojan  point 
orbit.  In  each  c&sa   the  approach  is  frost  the  negative  x 
direction,  and  she  vehicle  has  been  provided  with  approximately 
the  correct  relative  energy.  The  closest  point  of  approach 
to  the  Trojan  point:  is  considered  the  initial  point  for 
generating  thu  orbit:. 


* 


Fig.   '5.2       Effect  of  DispMcemsnt  Error  in  Establishing  a 
Tiojan  Point  Orbit 
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Under  these  conditions*  a  miss  to  the  left  of  the 
Trojan  point  results  in  an  orbit  whose  maximum  displacement 
from  the  Trojan  point  is  roach  smaller  than  that  of  the  orbit 
resulting  from  a  miss  to  the  right  of  the  Trojan  point.  In 
an  actual  attempted  injection  of  a  vehicle  into  this  type  of 
orbit  it  appears  we  should  evaluate  the  expected  angular 
error  in  guidance  .  and  then  aim  to  the  side  which  will  insure 
a  resulting  orbit  of  the  first  type. 
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CCETORISOH  OF  THE  ACCUAt  EQUATIONS  OF  MOTICS? 

^jltfinnn^jXMMuCXtBarjMtfw^n  wmuinwwa«»«tiB»»flW*<W!  I   ir    i"ln'i       i"i-i     T'w    i    -i--ii"ir       ■r~ri      i  niiMi  w  n  tt>  [fiTwr  rTri«~T 

WITH  THE  LI&&RI2ED  SQUATIC^S 


* 


We  will  define,  1S   /he  residual  force,  to  be  the 
vector  difference  between  the  actual  feree  exerted  oa  & 
vehicle  by  the  «oiabiised  gr*  vitational  and  centrifugal  force, 
and  the  fictitious  force  wkich  vould  be  exerted  en  the 
vehicle  by  the  linearised  licdel  of  the  system. 


■2* 


1c    *  g*ad  u  -  grad  u'  »  grad  u* 


wh^re 


u  »  u«u 


and  the  quantities  u  aid  u'  are  defined  by  equations  4,5 
and  4.6. 


2«£f«»i 


>  X 


Figure  6.1.     Effect  of  the  Residual  F^rce  on  the  Acceleratics* 
ols  tlvi  Vehicle 
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The  sre/ddual  force  directly  represents  the  effect  of 
all  higher  order  terms  of  the  equations  of  motion  which 
were  negl acted  to  their  linear isae ion.  As  snown.  in  Pig. 6.1, 
if  the  vehicle  U  assumed  to  follow  the  solution  of  the 
linearised  equations  ao  a  first  approximation,  then  the  resi- 
dual force  vector  show©  the  tendency  of  the  actual  sysfcda 
to  cause  the  vehicle  to  deviate  from  this  path* 

It  would  be  convenient  if  the  residual  force  field 
were  fcimd  to  have  some  significant  pattern.  For  example, 
if  a  solution  of  the  linearised  equations  was  formed  $?h@re 
the  residual  force  was  *i$s&y&   Inward  eowartutfte  Trojan  point, 
this  would  he  an  indication  of  the  existence  of  a  stable 
orbit  somewhere  within  this  region. 

Several  points  of  the  residual  Jaccbiaa  function*  u  , 

were  calculated  by  the  saethods  shown  £&  Appendix  8*  The 

results  of  these  calculations  are  shown  in  Fig.  s.l.  The 

residual  force  is  the  gradient  of  the  function  u",  i.e., 

it  acts  along  the  direction  of  greatest  rate  of  increase  of 

v'f  . _ 

u  with  distance,  and  its  magnitude  is  equal  to  that  yatft  of 

increase.  The  ploc;  of  Fig.  B.I  is  neither  sufficiently 
complete  nor  accurate  enough  for  a  numerical  determination 
of  the  residual  force.  In  my  case,  it  is  more  easily  calcu- 
lated diractly  by   component  &  with  £he  s&theds  of  Appendix  8. 
Instead  the  pic   i  meant  to  give  a  rough  indication  of  the 
direction  and  magnitude  of  the  residual  force  as  a  function 
of  displacement  from  the  Trojan  point. 
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No  significant  pattern  of  residual  force  distribution 
was  found.  If  we  trace  out  an  orbit  about  the  Trojan  point, 
we  find  that  the  residual  force  f lunctuates  in  magnitude 
and  direction  with  no  readily  observable  net  result.  We 
also  see  that  the  closer  we  are  to  the  Trojan  point,  the 
smaller  is  the  residual  force. 

So  far,  we  have  found  nothing  to  indicate  that  our 
linearised  model  is  not  a  good  approximation  of  the  actual 
system.  The  deviations  show  no  significant  patterns,  and 
they  seem  to  decrease  regularly  as  we  get  closer  to  the 
point  about  which  the  equations  were  linearized.  However, 
we  are  still  not  guaranteed  that  an  orbit  about  the  Trojan 
point  is  absolutely  stable  for  "small  disturbances",  nor  do 
we  have  any  information  about  how  small  the  disturbances 
must  be  to  obtain  "partial  stability,"  i.e.,  keeping  a 
vehicle  within  specified  limits  of  displacement  from  the 
Trojan  point  for  a  specified  period  of  time. 
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the  momm  of  vEsmmma  sxmu,vsn  by  digital 
cmmmsim  eg  the  orbit 

Usually  when  we  seek  an  analytical  solution  of  a  system 
of  non- linear  differential  equations?.,  we  feel  ghat  wa  can 
always  rail  back  on  numerical  method  is  all  else  fails.  We 
would  like  to  have  a  solution  in  closed  form  or.,  in  this  case,, 
an  analytical  upper  bound  on  the  orbit,  but,  in  general,  this 
is  a  convenience  which  saves  us  the  tim®  and  expense  of  resorting 
to  a  digital  computer*  Here,  there  is  seme  doubt  about  whether 
we  can  obtain  a  satisfactory  answer  to  the  stability  question 
by  numerical  methods* 

The  obvious  method  is  merely  to  integrate  the  equation 
of  motion  numerically  over  as  long  a  time  interval  as  is  required, 
and  thus  obtain  a  numerical  representation  o£  the  position  o£ 
the  vehicle  throughout  the  specified  time  interval.  The  limita- 
tions on  this  method  fall  into  two  classes* 

The  first  limitation  is  our  inaccurate  knowledge  of  the 
physical  constants  of  the  system.  Although  it  la  believed  that 
our  estimates  of  the  earth-moon  i    ratio  may  be  an  error  by 
as  much  as  0*3  per  cent,  our  present  linar  rocket  experiments 
should  significantly  advance  ©ur  knowledge  in  this  area.  In 
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any  case,   this  is  not  a  basic  limitation  of  the  invest!; 

of  stability*  If  we  can  obtain  She  answers  to  she  stability 

quest ions  based  on  a  model  of  the  ear£h«moon  system  where  -«^ 

is  exactly  0*01213,  we  will  have  taken  a  long  step  toward 
obtaining  answers  for  the  actual  system*  If  later ,  with  more 
precise  knowledge  of  the  mass  ratio,  we  feel  that  the  results 
do  not  apply  to  &he  actual  system,  the  procedures  for  compu- 
tation will  still  apply. 

The  major  limitation  on  determining  a  precision  orbit  is 
the  errors  inherent  in  the  method  of  computation.  The  effect 
of  the  choice  of  a  particular  method  of  computation  is  dis- 
cussed by  Baker  and  Others  <6}*  Basically,  these  errors  may 
be  classified  as  truncation  error  and  rounding  error.  Trunca- 
tion error  results  from  the  use  of  a  finite  interval  over  which 
the  value  of  the  integrand  is  approximated*  It  can  be  reduced 
by  decreasing  the  else  of  the  Interval,  and  thus  arsore  closely 
approximating  the  value  of  the  integrand  over  the  range  ot 
integration*  Rounding  error  is  &u®  to  the  limitation'  of  compu- 
tation to  a  finite  number  of  digits,  and  the  error  in  each  step 
of  the  ecmpu£s£iGB  is  propagated  throughout  the  entire  problem* 
The  resulting  rounding  error  then  increases  with  the  number  of 


These  errors  can  be  reduced  by  intelligent  selection  of 
methods  of  computation.  For  example,  for  m  orbit  of  this  type* 
it  is  felt  theft  cossputation  of  perturbations  from  a  reference 
orbit  would  be  more  efficient  than  dlrsct  Integration  of  the 
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equations..  However*  the  combination  of  truncation  error  and 
rounding  error  still  can  not  ba  eliminated*  For  any  precision 
orbit  computation*  there  is  a  practical  limit  to  she  duration 
of  the  orbit  which  we  can  compute  without  a  loss  of  validity 
from  the  etsnulatlve  error,  Thus^  for  a  general  lunar  probe, 
regardless  of  refinements  in  the  method  of  computation,  w& 
can  specify  where  it  will  be  on  its  first  pass  near  the  moon, 
we  can  estimate  where  it  will  be  a  few  months  later,  but  we 
have  no  idea  where  it  will  be  several  years  later. 

Intuitively,  it  is  felt  that  we  can  do  better  than  this 
with  a  periodic  orbit  such  as  ours.  It  seems  that,  if  the 
orbit  is  periodic,  we  can  compute  a  precision  orbit  to  complex 
one  revolution  about   the  Trojan  point,  and  the  orbit  is  stable 
if  it  passes  closer  to  the  Trojer*  point  the  second  tim®  around. 
But  we  heva  neglected  the  taction  of  the  vehicle  In  two  modes 
with  incommensurable  periods  ♦  A  hypothetical  vehicle  moving 
in  the  linearised  system  will  or  will  not  meet  the  above  cri- 
terion according  to  its  initial  conditions,  but  wa  know  that 
the  orbit  of  this  vehicle  is  stable.  Xfc  resialna  within  upper 
and  lower  bounds  determined  by  the  initial  conditions,  and  wil 
pass  arbitrarily  close  to  any  point  in  this  region  if  we  wait 
long  enough* 

The  nerst  thought  is  that  we  should  then  compare  the  computed 
precision  orbit  with  the  solution  of  the  linearised  equations • 
But  again  this  comparison  gives  no  conclusive  evidence  of  sta- 
bility* At  any  comparison  poiafe,  the  motion  of  the  vehicle  in 
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the  actual  system  will  fesvo  altered  she  proportions  of  rela- 
tive energy  between  the  two  modes.  We  may  find  that  Che  ac- 
tual orbit;  is  closer  to  the  Trojan  point,  but  if  we  compute  a 
new  linearized  orbit  using  an  initial  conditions  the  actual 
displacement  and  velocity  at  the  comparison  point,  we  may  find 
that  the  predicted  upper  bound  of  subsequent  displacement  fron 
the  Trojan  point  is  greater  than  that  predicted  by  the  first 
linearized  orbit* 

We  can,  of  course  calculate  the  predicted  upper  bound  and 
continue*  watching  this  upper  bound  for  an  increase  or  deerc- 
But  now  ve  are  back  to  the  problem  of  not  being  able  to  watch 
it  long  enough,  because  of  the  cumulattvs  error  of  the  numerical 
calculations*  Xt  is  this  cumulative  error  that  may  mask  the 
stability  effect  that  we  wish  to  determine,  especially  if  the 
instability  of  the  orbit  is  small •  He  am  working  in  a  regl 
of  cis lunar  space  roughly  bounded  by  two  ellipses  defined  by 
the  solution  of  the  linearised  equations*  Since  we  see  no  sig- 
nificant pattern  in  the  distribution  of  the  residual  forces  (or 
non- linear  effects)  In  this  region,  we  should  cover  the  entire 
region  rather  uniformly  in  our  calculations.  Otherwise,  we 
cannot  be  sure  that  the  nerct  few  revolutions  of  the  vehicle 
about  the  Trojan  point  will  not  pass  through  points  where  the 
non- linear  effects  are  entire!    ferent* 

The  length  of  &i*ae  necessary  for  an  orbit   to  pass  within 
a  specified  distance, o  ,   from  every  point  in  its  permitted  region 
varies  with  the  distribution  of  the  motion  between  the  two  modes 
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of  periodic  saotlon*  If  the  motion  of  the  vehicle  ia  priraari 
confined  so  one  sods  with  the  other  of  amplitude  no  greater 
than  6  ,  then  the  necessary  time  is  no  longer  than  the  period., 
of  the  primary  mode*  If  the  amplitude  of  the  secondary  mode 
is  large,  the  tints  may  be  prohibitive  due  to  cumulative  error* 

This  leads  us  to  consider  investigation  of  the  special 
cases  where*  based  on  the  linearized  model  of  the  system,  the 
initial  conditions  are  such  as  to  excite  only  one  mode  of 
vibration*  The  trouble  with  this  is  that  we  should  aspect 
the  non- linearities  of  the  actual  system  to  cause  transfer 
of  energy  to  the  other  mode*  with  a  cussulative  effect  until  we 
are  back  with  the  general  case  again* 

Xt  seems  possible  that,  we  may  obtain  useful  information 
hy  constraining  the  orbit  to  the  simple  ellipse,  and  inte- 
grating the  effect  of  the  residual  force  over  one  revolution 
about  the  Trojan  point.  Intuitively*  it  is  felt  that  a  net 
inward  force  would  be  an  indication  of  stability  while  a  net 
outward  force  would  indicate  instability.  If  this  were  done 
for  orbits  of  various  siaes*  we  would  then  h&ve  an  indication 
of  stability  versus  average  distance  from  the  Trojan  point* 
this  investigation  has  disclosed  no  mathematical  Justlflcati 
for  asserting  the  raeanisgfulness  of  this  "stability  indication". 
However*  It  is  felt  that  the  discovery  of  a  valid  stability 
criterion  will  require  calculations  of  this  &yp®* 
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CHAPTER  8 

LINEARIZATION  OF  XRE  EQUATIONS  OF  MOTION 
AT  POINTS  DISPLACED  FROM  THE  TROJAN  POINT 

8.1  Linearization  of  the  Equations  of  Motion  at  the  Position 
of  the  Vehicle 

In  a  digital  computation  of  the  orbit  of  a  vehicle  about 
a  Trojan  point,  it  seems  most  efficient  to  use  a  reference 
orbit  and  compute  deviations  from  this  orbit.  For  a  reference 
orbit,  the  first  choice  that  comes  to  mind  is  the  solution  of 
the  linearized  equations  as  developed  in  Chapter  5.  However, 
for  a  small  segment  of  the  orbit,  we  should  get  a  better  fit- 
ting reference  orbit  by  linearizing  the  equations  of  motion 
at  a  point  on  the  orbit,  rather  than  at  the  Trojan  point. 

Linearizing  equations  4.3  and  4.4  at  the  point  (a,b) 
rather  than  at  the  point  (0,0),  the  Taylor  expansions  are 

&-  »  *S  (a,b)  +  ^  (a,b)(x-a)  +  §-*   (a,b)(y~b) 
dx     dx        dx  dxdy 

(8.1) 


£H  .  *H  (a,b)  +  2Ji  (a,b)(x-a)  +  ^S(a,b)(y-b) 

dy   &y        dx^y  hy 
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For  convenience  in  writing  the  equations,  let  R*  and  R^ 
be  the  distances  from  the  point  (a,b)  to  the  earth  and  moon 
respectively*  Then  from  Figure  4.1 

R?  -  (b  + -)2  +  (a  *  k)2"  b2  +*3b  +  a2  +  a+l    (8.3) 

1  2 

JT  2  2 

R2  -  (b  -K3)  4-  (a  -  h)   •  b2  +  €b  *  a2  -  a  +  1    (8.4) 

2  2 

Equations  8.1  and  8.2  can  be  written  as 

$&  »  k,  +  k2  (x-a)  +  k3  (y-b)  (8.5) 

dx 


&  -  kA  +  k~  (x-a)  *  k,  (y-b)  (8.6) 

ax    '    J         J 


where,  from  equations  4.17  through  4.19,  with  substitution 
of  equations  8.3  and  8.4 

k,  =—  (a,b)  -  -e+aH  -A  (ah%)  R;3-€(a-%)  Rl3       (8.7) 
1   dx  A         * 

k2  =^-|  (a,b)  =  1-/*R~3-€R  +  3^(af?3)  R-,  *  3e(a-*>)  Ro 
dx 

(8.8) 

k  =lji  <a,fo)  =  3A(a**|)(b*-|)  R,  4-  3€(a-^)(b+  ~)  Ro 
J   dxdy  2   1  2   2 

(3.9) 
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^3 


/3x  ^ 


•    ^  (a,b)  mb  +  '&  -%(b  +  *1)  R.     -«  (b  *  ~)  rI     (8.10) 


/*  „'3 


.     2 


S  -  *8  ,..,,,  . 


-3 


« 


1-dR^  -  «B2"  ♦  3A(b  ♦  ~)  R^  +  3«(b  ♦  ■§)  R^ 


(3.11} 


Xn  llne&rlsed  fozsa*  eqts@fel.ons  4*3  sisd  4.4  can  be  bitten 


£    •  2  y  o  (k|  •  akg  -  fek^)  *  kg»  *  &3y  (8.12) 

y*  +  U  »  (k^  -  ek^  -  bk^)  *  k~K  *  &5y  (8.13) 

These  equations  can  be  wgt&tm  In  differential  operator  form 
as 

(p*  -  kg)  x  -  <2p  *  k^)  ywk^-  akg  -  bk3 
(2p  -  k3)  x  +  <p  -  k,-)  y  «  k^  -  akg  -  bk^ 


Hie  dbar&cteristle  eq^afcloia  of  this  system  Is 

<pa  ~  k2)  (pZ   -  kj)  *  (Sp  *  k3)  (2p  -  R3)  •  0 

ox* 


p^  +  (4  m   ^  „  &  )  p2  «$.  (kg%  -  kg)  «  0 
end  the  roots  of  the  characteristic  equation  ere  given  by 

r 


(k2k5  -kp 


so 


(t      II 

A  stable  solution  is  obtained  if 

<4-k2-k5)2  2 
J-1-  -k2k5*k^o 

Substituting  for  the  k's  from  equations  8.7  through  8.11, 
the  stability  criterion  is 

9  /  H    M    €   \2     i,  A    >     €,         9€/i(b  + .i) 

—  \  ^f    ***    "Try   "  2(  **qr  -f  — — *;   -  2 

4       R*  Rf  R^         r3  ^r5r5  7Q 

(8.14) 

Equations  8.12  and  8.13  can  be  solved  to  obtain  a  solu- 
tion  in  the  form  of  constants  plus  pure  sinusoids  if  the 
criterion  of  Eq»  8.14  is  met.  If  the  criterion  is  not  met, 
the  solution  contains  exponentially  increasing  sinusoids. 

Whether  or  not  any  significance  can  be  attached  to  this 
criterion  is  not  immediately  apparent:.  If  the  solution  resulting 
from  a  linearization  of  the  equations  about  a  point  displaced 
from  the  Trojan  point  is  to  be  considered  better  than  the  solu- 
tion of  Chapter  5,  it  can  be  used  only  over  a  small  segment  of 
the  orbits  and  then  a  linearization  in  the  vicinity  of  a  new 
point  must  be  accomplished.  We  cannot  then  categorically 
state  that  the  stability  of  a  specific  periodic  solution  is 
significant  if  we  consider  only  a  small  fraction  of  one  period 
of  the  solution. 

On  the  other  hand,  it  is  felt  that,  if  the  vehicle  passes 
through  a  region  where  the  solution  of  its  linearized  equations 
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is  unstable *  share  should  he  a  tendency  toward  stability  of 
the  complete  orbit.  Xt  seems  that  each  time  the  vehicle 
passes  through  this  region,  there  would  be  a  cumulative  ten- 
dency for  it  to  be  pushed  further  away  from  the  Trojan  point. 

Time  was  not  available  for  a  survey  of  the  region  in  the 
vicinity  of  the  Trojan  point  to  determine  where  the  criterion 
of  Eq.  8.14  is  met-  Roughly  it  seems  that  it  will  not  be  met 
for  either  x  or  y  greater  than  about  0.1.  In  any  case,  we 
should  first  have  a  better  justification  of  the  significance 
of  such  a  survey •  Although  a  time  limit  on  the  pressnt  study 
precluded  further  investigation  along  these  lines,  it  is  felt 
that  there  may  be  something  useful  in  this  area. 
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QgHSR  H&KTOS  OF  ORBITAL  MEMCOSm 

At  this  point,  the  ignorance  of  the  author  of  the  intri- 
cacies of  classical  celestial  mechanics  is  probably  already 
apparent-  The  literature  has  been  searched,  but  not  thoroughly 
digested*  1£  the  method  by  which  our  problem  will  be  solved 
remains  hidden  in  thi&  field,  it  oan  most  probably  be  found  in 
the  work  of  Aurel  Vintner  <7) . 

The  recent  Russian  vosU  on  the  subject  has  not  bean  tho- 
roughly examined  by  the  author ,  partly  because  of  the  non-aval 
bi'lity  of  trenslaticns*  The  work  of  Merman  (8)  is  bri&fly  dis- 
cussed  by  Leimanls  (9) •  This  method  seems  to  be  also  described 
by  an  abstract  of  a  paper  by  Sehigolev  (10); 

"The  motion  of  an  asteroid  is  calculated  according  to 

approximation  formulas :ie  approximation  values  desired 

axe  obtained  by  averaging.  The  average  siethod  of  the  author, 

however,  differs  from  the  well-known  method  of  celestial 
mechanics  inasmuch  as  the  author  does  not  average  the  power 
function,  but  certain  skillfully  chosen  parts  of  the  different 
tial  equations  and  thus  obtains  solvable  differential  equations". 

This  method  was  investigated  far  enough  to  ascertain  fchat 
periodic  solutioas  were  indeed  obtained,  but  not  far  enough  to 
judge  the  validity  of  these  solutions. 
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<mma  10 
cmcwB  ecus 

As  stated  la  Chapter  1,  *?©  do  not  have  specific  answers 
to  the  problem  under  consideration.  If  w©  are  askad  "Stow 
can  the  orbital  stability  of  a  proposed  teojan  poise  satellite 
be  predicted?'5,  we  csa  suggest  tha  following  approach* 

1.  E    ite  the  expected  guidance  inaccuracy  1st  terms 

of  velocity  and  d  .   ■■  at  the  T*    polar,* 

the  aisa&ag  point  as  sug&  './action  5*6  to  a&ti&il&e  the 
else  of  the  resulting  orbit  about  the  1!  i  pointy  based  ©a 
the  solution  of  the  linearised  equations » 

Z.     Compute  the  orbit  around  the  fro j  as*  poll    »lag  as 
incremental  reference  orblte,  solutions  of  the  equations  of 
Chapter  8*  Ths  length  of  the  incrjasants  used  la  the  eesipu* 
nation  will  be  picked  to  minimise  the  combination  of  truncati© 

error  «nd  rounding  error*  Also  m,  attest  should  be  sa&de  to 
predict  how  the  ®tm  of  this  error  will  is    *e  with  time. 
If  Che  computed  orbit  diverges  from  the  Troj&n  p©iM&  before  tfa 
©owputafcicm  error  become  stion  is  answer 

'    ■.■:■:'■■>■  '"-v-  qi:w&%   is  much  easier  to  obtain  Chan 
a  correct  poatlv-:     r.  If  t3  kit  diverges  fri- 

th© Trojan  point 


M 


this  is  due  Co  en  actual  Instability  of  the  orbit*  or  a  result 
©>:  the  accumulated  eosaputafcion  error* 

3*  Further  work  on  the  orbital  stability  of  Yrejan 
)  satellites  should  probably  be  directed  toward  the  problem  o£ 
identifying  and  reducing  the  effect;  of  eoiaputatioaal  errors* 

Also  it  is  felt  that  further  investigation  should  be  ctade 
of  the  significance  of  the  region  in  which  stable  solutions 
of  the  linearized  equations  are  obtained*  as  discussed  in 
Chapter  8* 
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APFEHDIX  A 


HUEIERlCftL  SCJLuTXQgg  01'  THE  LINEARIZED  STATIC;  IS  &   MQEXQI 


The  following  c&lculationif  were  made  to  permit,  calcu- 
lation of  a  solution  of  the  linearised  equations  dirasc.ly 
from  specified  initial  conditions.  For  numerical  worl^  a 
value  of  €  of  exactly  O.OJL213  i?as  used.  Althrugh,  the 
physical  quantity  in  the  earth  -moon  system  is  only  knoim. 
to  within  about  0.3  per  cent*  the  calcul//  n  were  car?; iec 
to  ten  significant  figures  to  permit  ui  s   J  tba  results  in 
subsequent  numerical  integration  comput; itions,  if  desired. 
The  results  are  then  consistent  with  .1  linearised  model 
of  the  system  with  the  mas;  ratio  specified  by  €  *»  0.01213 


/I    *  i  ^  i  »  0,98787 


l-27e  (a-  e )  0 ,32?.4735232 


pf  »  %  (1  ►  0.8224735232)  »  0,9112367616 


>6 


pf  -  h   (1-0.8224735232)  -  0,08-87632384 


*\ 


*>2 


»  0.9543872205 


-  0.2979316002 


1^  <X~2«)  =  1.267523442 


Substituting  inifcisX  conditions  x  (0),  y  (0).  &  (0) 
and  £  <0)  inc©  equations  5.1  and  5.2 

x  (0)  -  Bx  *  Dx  <A,1) 

y  (0)  »  Ba  +  D2  (A.2) 

*  <°)  °  *>1AX  *  *2C1  <A*3) 

y   (0)  *  pxA2  *  p2C2  (A.4) 

Combining  equations  5.3  and  5*4  we  can  get  A,  in 

terms  o£  B7  and  B« 
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apjAx 


tH> 


(1-2*)  tx  4   <pf  *    |)  B2 


<A«5) 


Similarly  fir/a  s^iafiioa©  5*5  and  5.6 


*2C3 


at      >  r mi    sx* 


Bl  *  <4  *  J>  D2 


Fran  e^fcioas  A»l  «£*!  A*2 


<A,6) 


%  ° 


(A.  7) 


D.j.  ®  ^  (5)  -  Bg 


(4.8) 


;kms  4*7  ©sa£  A.8  iskto  m*  4.6 


*«%  ar  np~  W*2* )  [»  <°)  -  h|  ♦  <P§  *  ~)  {jW-bJ 


«Bd  A*9 


.  7.*3  by  a 


<A.9} 
fsoa  squatlosifi  A*S 


(4,10) 
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2x  (0)  •    ^-^*    <l-2«)  x 
4 


B, 


2    2 
Pi  P2 


x{0)- 


p2*    £ 
r2        4 


<A.U} 


t 


<A.12) 


Subseifcafci&g  sususglc&l  values 

B2  *  -1*541111544  x<0)   -2*84357287?  y(0)  *  ?. 431(589219  x  (0) 


©£  A*,,,  «*  can  get 


BB 


3.3  32&I  5*4  &©  alimin&fce  A^  iastead 
2b  gesnis  of  B^  mz&  Bg 

3*3 


<1~2«)  B2     (A.  14) 


Similarly  fseea  ©quafciesaus  53  s&d  5*6 


i>/f 


r-     a      f» 


?  *  5)  8     ~    .23  <X~2s)  i>     (A.15) 

*      4      *  4  * 
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Substituting  equations  A.'/  and  A.U  into  Eq,   A.  15 

2?2%  •  -  (  l\  +  *)  [*<u)  ..  bJ-  Y  ^-2S>&<u)  -  B2j 


(A.  16) 

Multiplying  Eq.  A. 4  by  2  and  substituting  from  equations 
A.  14  and  A.  16 

2y  (0)  -  -  (pf  +  ^)  Bx  -  "  (l-2e)  |B2  +  y(0)  -  B2  J 

.(p|  +  |)    ^(0)  -  B  J 

(A. 17) 

2y  (.)  -  -  l^i  (1-26)  y(o)  -  (pf-p|)  Bx  -  (p§  *  |)  *(0) 

(A.  13) 

2   3  3>G_ 

bi  -  -  -^ttV  x(0)  -  =FI¥iI  y<u)-  -rV~  y<u) 

(A.  19) 

Substituting  numerical  values 

Bt  •  -1.019805762  x(0)  -1.541111544  y(o)  »  2.431689219  y(u) 

(A.20) 
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A^  caa  now  fee  found  from  uS 


(A. 21) 


Ax  -  0.6639118009  B^  +  1,653013473  32       (A.22) 

Substituting  from  equations  A* 20  and  A.  13 

Ax  *  3,228854338  x(0)  -  5-731588091  y(0) 

•14.026423771  i(0)  -  1*614427169  y<0j 

(A.23) 

Similarly  from  Eq.  A.  14 

2.3         3^3* 
L  +  4         ~5 — 

Tp^ —  Bl  "  — "^ 

(A.  24) 

A2  «  -0.8701335645  Bx  -  0.6639118009  B2 

(A.  25) 

Substituting  from  equations  A. 20  and  A.  13 


f^p.     tS     •»      i  m  "  ■  '   r  ii  i  I  Ot      **     i   i  m  .m'  ii  V  i.  i  <m  i  ufa 
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^  •  1*910529363  x(0)  *  3.228034338  y(Q) 

-1,614427169  x   (0)  +  2.115894408  y  (0}    <A.26) 

Substituting  Eq.  A, 20  into  Bq.  4*7 

t^  •  2.019805762  x<0)  *  1.541U1544    y<0)  +  2,431689219  y  (0) 

(A.27) 

Substituting  Bq.  A.  13  into  Bq,  A.8 

D2  *  1.541111544  x  (0)  +  3.843572677  y  (0)  -  2.431689219  k  (0) 


C,  can  be  found  fjrata  Eq*  A,  6 

lp  (1.2t)  p|  *    | 

■Ml  - 

2p 


1  2o.  -  ^2  2 


2 


(A.29) 
Cx  o  2.127205440  T*x  -¥  3,925000296  I>2  (A, 30) 
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Substituting  from  equations  A. 27  and  A.2U 

Cx  -  10.3454050/  x<0)  *  Li. 364234/6  y(o) 

-9.544330904  s(0)  +  5.172702535  y(o) 

(A. 31) 

Similar3.y  from  Eq«  A.  15 

2   ,      3 
C2  =  "  2p/  Dl" 

(A.32) 

C2  «  -  1.407643965  Dx  -  2.127205440  D2  (A. 33) 


Substituting  from  equations  A. 27  and  A. 23 

C2  -  -  6.121420251  x(0)  -  lo. 34540507  y(o) 

*  5.1/2702535  x  <o)  -  3.422952654  y(u) 


(A. 34) 
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Given  specified  initial  conditions,,  the  constants 
of  integration  casa  301?  h®  calculated  frca  equ&tlcna  A.X3-, 
A*20#  A*23,  A- 25,  A*2?5  A*28,  A»31  and  A.34*  These  equa- 
tions wars  cross-checked  numerically  hy  other  formulas  end 
found  to  be  consistent* 

Xh©  substitution  of  number©  into  this  equation  tells 
us  little  about  what  sort  of  an  orbit  we  would  get  with 
slightly  different  Initial  /.©editions,  as  does  the  raathod 
of  Section  5*5.  Sostsver,  i'or  calculation  of  a  specific 
sasrbit.,  the  n$finarlea,3.  wor&  ia  ssuch  simpler*  If  this  s©Xu«» 
tioa  were  required  a  11  a  ralerenae  orbit  for  a  digital 
ccfsputation  of  the  orbit,  the  sbo*<?e  eftpatlofas  would  be  used 
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AggggggBC  B 

gJUMF/.ICAL  COMFi  TATXO^  OF  THE  RESIDUAL  JACOBIM  FPKCraOii 

To  «iet€7mtee  sobs©  msaerieal  values  of  u  axid  u*, 

equations  4»5  sad  4.20  *<ra«e  fcxa»s£©2med  to  pol(sx  coordi* 
t^ees  by  t*D  equations 


■  r  so& 


e 


y  *  p  sis  v? 


Tha  equations  fchcsa  becosae 

u  **  «.  epeoad*  %  (p2sln2d*^p^.»0+p2  coa20*oeos0+i) 

(sa) 

u«  «    i  +  1  p2co«i20+    ~"~  (3L*2«)p28ln  &eog6*«~  p2©3^  6 
2      8'  4  '  &  > 


^:~ 


Using  the  identities 


&+  cos2  & 


nist  v  *  co@    C7  »  l 


—  ©in  0  4-  %  cog  $«  gin  60°  sin  C^*-  cos  60°  cos  & 
2 


vT5 


co®  {#~6<r) 


2     sin  v-  %  eos  fcr  »  -  c©8  (  ft*  60°) 


the  equations:  £02?  u  and  n  •  become 


u  ■  -epcosd*  %  K>*  +  2j^cos  (  0-60°)  +  X  J 

+^{j?  *  2fco*  (8-60°)  *  XJ    %  e  W^  ^s0*K#*  60°}+; 


(B.3) 


U8CT    i  +    1   p2  +  Jl3  (i.2«)  P  2ain  2  0»  1  $2  cos  2  © 


(8.4) 


Usirag  a  vain©  foa?  e  of  0.01213,  the  emsntities  u  assd  u1 
calculated  for  0.3,  rjf      fxmi  0  to  0.5  «nd  fox 
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30°  iac?a&€S&@  o£  $  *    The  results  see  listed  in  Tab  Is  B*l» 
the  areei&ial  function  n"  «  u  -  u*  at  fchea©  points  is  shown 
in  Fig*  B#l« 
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Fig.  B.  1»  Approximate  Values  of  the  Residual  Jacobian  Function 
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TABLE  B.X 


TABULATION  flg  JACOBIAH  FUNCTIONS 


e 

f 

u 

u! 

e 

f 

14 

u! 

0° 

.1 

1.5041 

1.5033 

180° 

11 

1.5033 

1.5038 

-2 

1.5177 

1.5150 

»«£ 

1.5112  1.5150 

.3 

1.5429 

1.5333 

f3 

l  1.5803  1.5338 
1  1.5273  1.5600 

.4 

1.5325 

1.5600 

.4  ' 

•5 

1.6297 

1.5938 

.5 

1.5310 

1.5938 

30° 

•  1 

1.5108 

1.5111 

210Q 

.1 

1.5114 

1.5111 

.2 

1.5429 

1.5445 

.2 

1.5470 

1.5445 

.3 

1.5919 

1.6000 

|  .3 

1.6080 

1.6000 

.A 

1.6592 

1.6773 

.4 

1.6944 

1.6778 

.5 

1.7427 

1.7778 

.5 

1.8023 
1.5159 

1.7778 
1.5149 

60° 

.1 

1.5140 

1.5X49 

240° 

.1 

.2 

1.5529 

1.5595 

.2 

1.5692 

1.5595 

.3 

1.6134 

1.6338 

.3 

1.6717  1.6338 

.4 

1.6930 

1.7338 

.4 

1.8428  1.7378 

.5 

I  1.7927 

1.8716 

270° 

.5 

2.1177  1.8716 

90° 

.1 

1.5109 

1.5113 

.1 

1.6116 

1.5113 

.2 

1.5434 

1.5450 

1.5476 

1.5430 

•3 

1.5930 

1.6013 

.3 

1.6093 

1.6013 

•4 

1.6611 

1.6300 

.4 

1.6966 

1.6300 

•5 

1.7455 

1  1.7813 
1.5039 

•5  | 

1.8058 

1.7813 

120° 

.1 

1.5043 

300° 

*  X 

1.5035 

1.5039 

.2 

1.5183 

1.5156 

.2  | 

1.5119 

1.5156 

•3 

1.5441 

1.5350 

.3 

1.5219 

1.5350 

■ 

1.4 

1.5846 

1.5622 

.4 

1.5310 

1.5622 

•5 

1.0329 

1.5002 

1.5972 

.5 

1.5369 

1.5972 

150° 

.1 

1.5001 

330°  1  .1 

1.5002 

1.S001 

•2 

1.5010 

1.5006 

.2 

1.5012 

1.5006 

.3 

1,5049 

1.5012 

r  .3 

1.5041 

1.5012 

♦4 

1.5103 

1.5022 

}  .4 

1.5108 

1.5022 

■  ■--..u.- 

•5 

1.5222 

1.5034 

J  *5 

1.5229 

1.5034 
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